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ABSTRACT 

Multicarrier  (MC)  CDMA  systems  have  been  developed  and 
rapidly  gained  popularity  as  they  capitalize  on  both  MCM 
(multicarrier  modulation )fs  resilience  to  MUI  and  direct- 
sequence  ( DS )  CDMA's  robustness  against  frequency  selec¬ 
tivity  to  mitigate  both  MUI  and  the  inter-symbol  interfer¬ 
ence  (ISI)  caused  by  time  dispersive  channels.  However, 
the  well-known  near-far  problem  in  a  multiuser  setting  still 
places  fundamental  limitations  on  the  performance  of  MC- 
CDMA  communication  systems ,  and  this  issue  has  been 
considered  by  very  few  papers.  In  this  paper,  the  theoretical 
near-far  resistance  of  the  MMSE  detector  in  the  uplink  MC- 
CDMA  (with  and  without  cyclic  prefix  (CP))  setting  is  de¬ 
rived,  It  turns  out  that  the  near-far  resistance  ofMC-CDMA 
without  CP  has  the  same  form  as  that  of  the  DS-CDMA  ex¬ 
cept  that  the  user  codes  are  IFFT  transformed.  Our  formu¬ 
lation  for  the  near-far  resistance  is  applied  to  either  blind  or 
non-bljnd  MMSE  detectors .  Computer  simulations  confirm 
theoretical  findings, 

1.  INTRODUCTION 

New  multimedia  applications  in  the  field  of  mobile  com¬ 
munications  require  high  data  rate  transmission  under  the 
frequency  selective  fading  environment.  The  current  DS- 
CDMA  systems  may  become  impractical  when  the  data  rate 
reaches  very  high  speeds.  The  very  short  chip  duration  leads 
to  severe  interchip  interference  (ICI)  due  to  multipath  fad¬ 
ing  channels  and  difficulties  of  synchronization.  The  MC- 
CDMA  technique  does  not  suffer  from  such  shortcomings. 
As  a  hybrid  of  multicarrier  modulation  (MCM)  and  DS- 
CDMA,  each  chip  of  the  desired  user’s  signature  sequence 
is  made  longer  and  is  modulated  onto  a  separate  subcar¬ 
rier,  which  has  narrow  bandwidth  undergoing  frequency-flat 
fading.  Because  synchronization  is  very  important  for  per¬ 
formance  of  MCM  approaches,  most  researchers  in  MC- 
CDMA  area  choose  to  focus  on  the  synchronous  trans¬ 
mission  scenario,  i.e.,  downlink  case.  However,  the  well- 
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known  near-far  problem  in  a  multiuser  setting  (in  uplink) 
still  places  fundamental  limitations  on  the  performance  of 
MC-CDMA  communication  systems,  and  this  issue  has 
been  considered  by  very  few  papers.  In  this  paper,  die  the¬ 
oretical  near-far  resistance  of  the  MMSE  detector  in  the  up¬ 
link  MC-CDMA  setting  with  or  without  cyclic  prefix  (CP) 
is  derived.  It  turns  out  that  the  near- far  resistance  of  MC- 
CDMA  without  CP  has  the  same  form  as  that  of  the  DS- 
CDMA  except  that  the  user  codes  are  IFFT  transformed. 
Our  formulation  for  the  near- far  resistance  is  applied  to  ei¬ 
ther  blind  or  non-blind  MMSE  detectors. 

The  notation  used  in  paper  follows  usual  convention: 
Vectors  are  denoted  by  symbols  in  boldface,  (•)*,  (*)T,  (.)H 
are  complex  conjugate,  transpose  and  conjugate  transpose 
of  (*),  respectively. 

2,  SYSTEM  MODEL 

Consider  an  asynchronous  MC-CDMA  system  with  N  sub¬ 
carriers  and  J  active  users  in  a  multipath  fading  chan¬ 
nel  Random  binary  (±1)  codes  of  length  N  are  used  as 
spreading  codes  to  distinguish  different  users.  We  will  be 
concerned  with  digital  wireless  transmissions  through  lin¬ 
ear  time  invariant  (LTI)  discrete-time  baseband  equivalent 
channels.  The  baseband  MC-CDMA  system  model  for  the 
jth  user  is  presented  in  Figure  1  on  the  next  page.  With  re¬ 
ceived  baseband  signal  sampled  at  chip  rate  1/TC  where  Tc 
is  the  chip  duration,  the  discrete-time  received  sequence  is 
given  by 

L 

fj(n)  =  ^2  hj(l)uj(n  - 1)  +  Vj (n)  (1) 

1=0 

where  fj(n),  %(n),  hj(n)  and  vj(n)  are  jth  user’s  re¬ 
ceived  signal  transmitted  OFDM  symbols*  chip  sampled 
multipath  channel  and  additive  white  Gaussian  noise,  re¬ 
spectively.  In  the  following,  we  will  only  consider  such 
discrete-time  equivalent  models  which  will  be  assumed  to 
satisfy  the  following. 

Assumption  1;  All  user  channels  are  FIR  of  order  <  L  - 
L  +  £>,  where  L  is  the  maximum  delay  spread  of  the  chip 


AWGN 

Fig.  1.  Baseband  MC-CDMA  system  model  for  the  jth  user 


sampled  multipath  channels  and  D  is  the  maximum  relative 
asynchronous  delay  among  the  users. 

The  MC-CDMA  model  with  CP  was  first  developed  by 
combining  the  OFDM  and  CDMA  techniques  [1](2].  How¬ 
ever,  due  to  the  CP  overhead,  researchers  are  more  inter¬ 
ested  in  detectors  for  the  MC-CDMA  models  without  CP 
[4]  [5].  We  will  first  derive  the  discrete-time  MC-CDMA 
model  without  CP  Group  the  serial  %(n)  into  blocks  of 
N  (N  »  L)  (with  CP  the  block  size  is  JV  +  L).  De¬ 
fine  the  jth  user’s  ith  transmitted  block  to  be  Uj(i)  = 
[v,j(iN).  uj(iN  +  1),  •  •  ■ ,  Uj(iN  +  N  -  l)f  and  f  j(i)  = 
[fj(iN),rj(iN  +  1),  •  •  • ,  fj(iN  +  N  —  1)]T  as  the  ith  re- 
ceived  block.  Using  (1),  we  obtain 


(i)  =  Hjo  Qj  (*)  +  Hji  uj  (i  -  1)  +  vj  ( i )  (2) 


where  Vj(i)  is  the  corresponding  jth  user’s  noise  vector, 
and 
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Note  that  depending  on  the  values  of  N  and  £,  there  may 
be  more  inter-block  interference  (IBI)  terms  and  corre¬ 
spondingly  more  5  lj  >  1.  In  the  MC-CDMA  sys- 
tem,  u  j(i)  =  WH  [cjU-  •  ■,c]N)t  bj(i)  =  WHc}bj(i)  = 
c'  j  bj(i)  where  WH  is  the  N  x  N  dimensional  IDFT  ma- 
trix  with  its  (m,n)ih  element  as  ^ 

Cj  and  bj(i )  are  the  spreading  code  vector  and  ith  trans¬ 
mitted  bit  of  the  jth  user,  respectively.  d  j  —  WH  c j  is  the 
IDFT  transformed  code  vector.  Define  the  (N 4- L)  x  (L  -h  1) 
matrix  Cfj  with  shifted  column  vector  c f  j  and  jth  user’s 


channel  vector  h;  as 


C'i  = 


C  j\ 


■  jN 


C'jl 


■  jN 


,hj 


MO) 
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After  carefully  examining  the  structure  of  Hji,  and  Uj  (i), 
we  have 


=  C j(ljN  +  1  :  IjN  +  N,:)  hj  bj(i) 

=  C'i(lj)h,6J(i)  =  Hi(/j)ft,(i)  (3) 

where  in  the  above  equation,  we  use  MATLAB  represen¬ 
tation  to  denote  C'j(lj)  as  the  submatrix  of  Cfj  from  row 
IjN  +  1  to  IjN  4-  N.  In  the  receiver,  the  received  signal  is 
first  feed  through  an  FFT  block.  The  total  received  signal  at 
the  receiver  is  the  superposition  of  the  data  signals  of  the  J 
users,  given  by 


j 

r(*)  =  EWM 

i=i 


where  W  is  the  FFT  matrix.  Since  the  receiver  does  not 
know  the  starting  time  of  each  desired  user  symbol,  then 
similar  to  a  chip  level  detector  [6]  for  an  asynchronous  D$- 
CDMA  system,  a  smoothing  factor  M  is  needed  to  capture 


a  complete  desired  user  symbol.  Let  Lh 
from  (2)-(4),  we  obtain 


=  maxL-  +1  and 
i= i  3 


XAf(n)  =  %b(n)  +  v(n) 

(5) 

where  Xa/W  =  [  rT(n)  rT(n  +  M-l)' 

T 

9 

b(n) 

=  hT  (n  -  Lh  +  l)  br(n  +  M-l) 

T 

b(n) 

=  [  bi(n)  •••  bj(n)  ) 

T 

v(n) 

=  [vT(n),...!vT(n  +  M-l)]T, 

'WH(L4-1)  •••  WH(0)  0 

1 

n  = 

0  ■■■  WH(ti  — 1)  ■■■  WH(0) 

and  H(Z)  =  [Hi(Z),  •  •  • ,  Hj(/)]t.  It  is  not  hard  to  find  that 
(5)  has  the  similar  form  as  that  of  the  DS-CDMA  [6],  Since 
the  received  signal  amplitude  is  the  combination  of  channel 
attenuation  and  the  transmitted  signal  amplitude  b(n),  we 
can  combine  these  two  effects  into  one  diagonal  amplitude 
matrix  A.  Thus  (5)  can  be  rewritten  as 

XM{n)  =  HAb(n)  +  v(n)  (6) 

with  b(n)  having  unit  norm,  and  the  columns  of  the  channel 
matrix  "H.  being  all  normalized. 

For  the  MC-CDMA  with  CP,  after  the  FFT  and  the  elim¬ 
ination  of  the  CP,  (4)  can  be  rewritten  as,  in  the  absence  of 
noise, 

ij{i)  =  W Rcp  Hjo  Tcp  Cj  bj (i) 

+  W  Rcp  Hj!  Tcp  WH  Cj  bj  (i) 

where  Tcp  and  Rcp  represent  the  processes  of  adding  and 
elimination  the  CP  [1],  Due  to  the  effect  of  Tcp  and  Rcp, 
Rcp  H ji  Tcp  =  0,  the  IBI  term  is  thus  eliminated.  Further¬ 
more,  Rcp  Hj0  Tcp  =  H;o  becomes  a  circulant  matrix  with 
first  column  as  jhj,  0  *  *  * ,  oj  . 

Fact  1  (Diagonalization  of  a  circulant  matrix)  [1];  An 
NxN  circulant  matrix  H^o  can  be  diagonalized  by  pre-  and 
post-multiplication  with  N-point  FFT  and  IFFT  matrices. 
Its  diagonal  elements  are  the  corresponding  FFT  transform 
of  the  first  column  of  the  circulant  matrix. 

Therefore,  the  resulting  channel  matrix  %  for  MC- 
CDMA  with  CP  is  reduced  to  H 

f  cnFi(l)  C/iFj(  1) 

n=  ■.  :  (7) 

cin  Fi(N)  cjpj  Fj{N) 

where  [Fj (1),  •  •  • ,  Fj (N)f  =  W  [hj,  0  •  •  • ,  o] T. 

The  following  assumptions  will  be  made  throughout  this 
paper. 

1.  The  symbols  bj(n)  are  uncorrelated  in  time;  6,(n)  and 
bj(n),  i  ^  j  are  uncorrelated. 

2.  The  channel  matrices  Ti  in  (6)  and  T-t  in  (7)  are  of 
full  column  rank  (known  as  the  identifiability  condition 
in  the  blind  multiuser  detection/equalization  literature 
[6][7]). 

3.  NEAR-FAR  RESISTANCE  OF  MC-CDMA 

In  this  section,  first  we  prove  that  fTOmse  =  ,4dR+Hd  is  the 
general  MMSE  detector  for  the  signal  model  (6),  where  “+” 
represents  pseudoinverse,  R  is  the  autocorrelation  matrix  of 
the  received  signal  vector  xm(h),  Hd  is  one  of  the  columns 


in  H  corresponding  to  the  desired  user’s  symbol  we  want 
to  detect,  and  Ad  is  the  received  signal  amplitude  of  the 
desired  user.  According  to  [10],  the  MMSE  detector  can  be 
obtained  by  the  Wiener-Hopf  equation 

f mmse  =  R  E  { Xm 

=  R-1WA[G,.-.!O,1)Oj--.,0]t 
=  AdR-1  Hrf  (8) 


where  bd(n)  is  the  desired  user’s  transmitted  symbol,  and  1 
is  in  the  dth  position.  Since  the  zero  forcing  (ZF)  detector 
will  be  used  in  the  derivation  of  the  near-far  resistance,  we 
first  show  that  when  noise  approaches  to  zero,  these  two 
detectors  are  proportional  to  each  other.  Assume  there  is  no 
noise  and  the  zero  forcing  detector  is  f zf,  Applying  to 
(6),  we  require 

f "fXM{n)  =  f^f"WAb(n)  =  [0  •  *  •  0  q  0  •  •  •  0]  Ab(n) 
=  aAdbd(n)  (9) 


where  a  =  f  ^Hd  is  a  constant  and  at  the  dth  posi¬ 
tion.  Then  right  multiply  (9)  by  x|(n)  and  take  expecta¬ 
tion,  we  obtain  fffE  [xAr(n)X^(n)j  =  aA2d Hf.  There¬ 
fore,  the  zero  forcing  detector  has  the  expression  f zf  = 
aA^R-1Hd.  It  is  straightforward  to  show  that  for  the 
noiseless  case  the  ZF  and  MMSE  detectors  are  related  by 
f  zf  ~  Ol  A.d  fmmse  • 

By  applying  the  zero  forcing  detector  to  the  received  sig¬ 
nal  vector  xm  (n),  the  output  contains  only  the  useful  signal 
and  ambient  Gaussian  noise.  The  energy  of  the  useful  sig¬ 
nal  at  the  output  is  Es  =  A2d{^fHdJlf  fzf,  the  variance  of 

the  noise  is  En  =  f  ff$zf  where  a2  is  the  power  spectral 
density  of  white  Gaussian  noise.  Using  the  definition  in  [8], 
the  asymptotic  multiuser  efficiency  (AME)  for  the  desired 
user  is 


ej(v)  **[Q-l(Pdl*))Y 

m  =  hm  =  hm  - ± -  - =L 

ar-+0  Aj  ff-40  Ag 


!  ff;  Ff  A\ 


f£m3gHdH  f  f  mmse 

,  7h  } 

*  mmse  1  mmse 


__  Hm  H f  (HA  g"  +  l) +  Hd  Hf  (HA  H»  +  l)  +  Hd 
c'-*o  Hf  (nAUH +<t21)+('HA'Hh +  <T21)+Hd 
=  Hf  ^w+A~1^+HdHf  HIt  +  A~1H+Hd 
Hf  nH+  A~1n+HH  +  A-^+Hd 

4— 2 

_ _ ‘Jjl _ 

Hf  n+H  A-1  A-1  W+Hd 

=  (10) 


where  Pd(o)  =  Q  •  an^  Q  *s  the  complemen¬ 

tary  Gaussian  cumulative  distribution  function.  In  (10)  we 


have  used  the  facts  that  (' HPMh)+  =  (H+)h  A ~lrH+, 

(' HhH)+  =  n+  (H+f  and  W+Hd  =  [0  -  -  010  -0]", 

where  1  is  in  the  dth  position.  From  (10),  it  is  seen  that  the 
AME  does  not  depend  on  the  interfering  signal  amplitudes. 
Thus,  it  is  equal  to  the  near- far  resistance  Vd  [8].  It  then 
follows  that 


%  =  Vd  = 


(id 


We  can  see  that  the  near-far  resistance  of  MMSE  detector 
of  MC-CDMA  without  CP  has  the  same  form  as  that  of 
DS-CDMA  [6]  [9],  The  only  difference  is  that  in  the  MC- 
CDMA  setting  the  user  codes  are  IFFT  transformed.  Due 
to  the  property  of  IFFT  matrix  (i.e.,  if  the  original  users’ 
codes  are  orthogonal,  the  resulting  IFFT  transformed  users’ 
codes  are  still  orthogonal),  from  (11),  one  can  expect  that 
the  near-far  resistance  of  MC-CDMA  without  CP  should  be 
quite  similar  as  that  of  DS-CDMA  systems. 

Following  the  steps  of  (10),  it  can  be  shown  that  the  near- 
far  resistance  of  MMSE  detector  of  MC-CDMA  with  CP 
has  the  similar  form  as  (1 1)  with  M  now  substituted  by  H. 
Since  the  multiuser  interference  (MUI)  is  reduced  through 
the  IBI  free  transmission  in  MC-CDMA  with  CP,  we  can 
expect  that  MC-CDMA  with  CP  is  more  near-far  resistant 
than  MC-CDMA  without  CP  at  the  expense  of  low  band¬ 
width  efficiency.  Finally,  we  should  point  out  that  (11)  is 
the  theoretical  value  of  near-far  resistance  and  may  not  be 
achieved  in  practical  situations. 


4.  SIMULATIONS 

Simulation  examples  are  presented  to  compare  the  theoret¬ 
ical  and  practical  near-far  resistance  of  two  blind  detection 
method  i.e,,  the  subspace  method  [6]  and  the  linear  predic¬ 
tion  method  [7]  under  both  DS-CDMA  and  MC-CDMA  set¬ 
tings.  The  practical  near-far  resistance  of  the  desired  user  is 

||fH  ||2 

calculated  as  r\  —  ■  fuf  in  the  noise  free  situation.  1000 
symbols  are  used  to  estimate  various  detectors.  In  all  of  the 
simulations,  the  users’  spreading  codes  were  randomly  gen¬ 
erated  binary  (dbl,  with  equal  probability)  sequences,  the 
number  of  subcarriers  N  =  64,  The  symbol  length  is  also 
64  samples.  The  multipath  channels  for  each  user  have  ten 
paths  which  may  not  be  uniformly  spaced  in  time.  The  max¬ 
imum  multipath  delay  spread  of  the  channel  is  one  fourth  of 
the  symbol  length.  The  path  gains  are  generated  according 
to  Gaussian  distribution  with  zero  means  and  unit  variances. 
The  maximum  asynchronous  user  delay  D  is  one  half  of  the 
symbol  interval.  In  DS-CDMA  the  processing  gain  is  cho¬ 
sen  as  64  to  retain  the  same  bandwidth  as  that  of  the  MC- 
CDMA,  The  channels  for  each  user  were  randomly  gener¬ 
ated  in  each  of  the  100  Monte  Carlo  runs.  All  input  symbols 
are  drawn  from  a  BPSK  constellation  and  then  multiplied 
by  various  magnitude  factors  to  generate  the  near-far  situ¬ 
ations,  All  users  other  than  the  desired  user  have  the  same 


Near-far  ratio 

OdB 

10  dB 

20  dB 

DS-T 

0.3127 

0.3108 

0.3073 

MC-T 

0.3255 

0.3272 

0,3089 

DS-P-SS 

0.2666 

0.2780 

0.2503 

MC-P-SS 

0,2693 

0.2821 

0.2411 

DS-P-LP 

0.2545 

0.2708 

0.2802 

MC-P-LP 

0.2501 

0,2853 

0.2886 

Table  1.  Comparison  of  theoretical  and  practical  near-far 
resistance  of  [6] [7]  under  different  near-far  ratio,  DS;  DS- 
CDMA,  MC:  MC-CDMA  without  CP,  T:  Theoretical  value, 
P;  Practical  (simulation)  value,  SS:  Subspace  method  [6], 
LP:  Linear  prediction  method  [7]. 


signal  power.  The  near- far  ratio  is  defined  as  20  log10 
for  i  f  1,  where  subscript  1  denotes  the  desired  user.  The 
smoothing  factor  M  —  2  is  used  for  all  algorithms.  For  the 
signal  subspace  in  [6],  the  rank  was  determined  by  taking 
empirically  selected  thresholds . 

Example  1:  MC-CDMA  without  CP  vs.  DS-CDMA 

In  this  example,  we  use  computer  simulations  to  verify 
our  conclusions  on  the  near-far  resistance  of  MC-CDMA 
without  CP  by  comparing  with  the  theoretical  values  (11) 
and  with  those  of  DS-CDMA,  The  simulation  is  set  under 
three  different  near-far  ratios  (i,e„  OdB  10  dB  and  20  dB). 
The  number  of  active  users  is  30.  The  results  are  presented 
in  Table  1.  From  the  table,  it  can  be  seen  that  the  near- far 
resistance  of  MC-CDMA  without  CP  is  almost  the  same  as 
that  of  DS-CDMA. 

Example  2:  Near-far  resistance  vs.  system  load 

In  this  example,  we  compare  the  practical  near-far  resis¬ 
tance  of  different  detectors  with  respect  to  different  number 
of  active  users  at  20  dB  near-far  ratio.  It  can  been  seen  from 
Figure  2  that  the  theoretical  as  well  as  the  practical  near- far 
resistance  of  DS-CDMA  and  MC-CDMA  without  CP  is  still 
quite  similar,  while  MC-CDMA  with  CP  has  much  better 
performance.1  As  mentioned  before,  this  performance  gain 
is  achieved  at  the  expense  of  reduced  bandwidth  efficiency. 
Furthermore,  it  can  be  seen  that  the  near-far  resistance  de¬ 
creases  when  the  system  load  increases.  This  is  because 
more  active  users  lead  to  more  MUI  and  thus  more  severe 
near- far  problem. 

Example  3:  Near  far  resistance  vs.  symbol  length 

In  this  example,  we  compare  the  practical  near-far  resis¬ 
tance  of  different  detectors  with  respect  to  different  length 
of  symbol  record  at  20  dB  near-far  ratio.  The  system  has  30 
active  users.  From  the  simulation  results,  it  can  be  seen  that 
when  more  number  of  symbols  were  used  to  estimate  the 

1  Linear  prediction  method  [7]  is  not  applicable  to  MC-CDMA  with  CP 


Fig.  2.  Near-far  resistance  vs  different  number  of  users, 
20  dB  near-far  ratio,  o;  Theoretical  near  -far  resistance  of 
DS-CDMA  (DS);  x ;  Linear  prediction  under  DS-CDMA 
(DSLP);  v:  Subspace  under  DS-CDMA  (DSSS);  □:  The¬ 
oretical  near  -far  resistance  of  MC-CDMA  without  CP 
(MNCP);  *:  Linear  prediction  under  MC-CDMA  without 
CP;  A;  Subspace  under  MC-CDMA  without  CP;  +:  The¬ 
oretical  near-far  resistance  of  MC-CDMA  with  CP  (MCP); 
<:  Subspace  under  MC-CDMA  with  CP  (MCPSS), 

detector,  the  performance  on  near-far  resistance  becomes 
slightly  better  in  all  simulated  algorithms  since  detector  f 
is  estimated  more  accurately; 

5.  CONCLUSIONS 

The  theoretical  near-far  resistance  of  MMSE  detector  in 
MC-CDMA  system  has  been  derived.  Numerical  compar¬ 
ison  of  near-far  resistance  between  DS-CDMA  and  MC- 
CDMA  is  presented.  It  is  shown  that  the  near-far  resis¬ 
tance  of  MC-CDMA  without  CP  is  very  similar  to  that  of 
DS-CDMA  systems,  while  MC-CDMA  with  CP  has  much 
better  performance  on  near-far  resistance  at  the  expense  of 
reduced  bandwidth  efficiency. 
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